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Abstract. We point out an apparently overlooked consequence of the boundary conditions 
obeyed by the electric displacement vector at air-metal interfaces: the continuity of the normal 
component combined with the quantum mechanical penetration of the electron gas in the air 
implies the existence of a surface on which the dielectric function vanishes. This, in turn, leads 
to an enhancement of the normal component of the total electric field. We study this effect 
for a planar metal surface, with the inhomogenous electron density accounted for by a Jellium 
model. We also illustrate the effect for equilateral triangular nanoislands via numerical solutions 
of the appropriate Maxwell equations, and show that the field enhancement is several orders of 
magnitude larger than what the conventional theory predicts. 
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1 INTRODUCTION 

Electromagnetic field enhancement in nanoplasmonic structures plays an important role in a 
number of applications such as nano antennas ||Tl|2l, single-photon emitters |3|, and surface- 
enhanced spectroscopy techniques in general I?], including surface-enhanced Raman spec- 
troscopy (SERS) JS]. Field localization is the key to understanding and further enhancing the 
resolution of SERS, and the use of highly engineered surface plasmonic structures has been 
envisioned |l6]|7l. Strong field enhancement has so far been inherently connected to singulari- 
ties m. 

It is commonplace to study the plasmonic response using Maxwell's equations with the bulk 
dielectric function of the metal, with an abrupt change at the air-metal interface 13 . Though 
clearly neglecting surface and confinement effects for the electrons in the metal, this is a cel- 
ebrated approximation which predicts the existence of surface plasmons forming the heart of 
nanoplasmonics. However, when samples involve true nano-scale feature-sizes, the validity of 
this approximation is challenged and new phenomena, such as non-local response and spatial 
dispersion ifTOl - fTTi and inhomogeneous electron densities due to quantum confinement 11811191 
are expected to play an important role. 
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In nanoscale metallic structures the electrons are subject to quantum confinement, which 
leads to a number of spectacular effects, such as the quantized conductance of atomic quantum- 
point contacts {20l , Friedel oscillations and confinement of electrons to quantum corrals on 
noble metal surfaces II2TI . and a size-dependent structure of nanocrystals 1221 . Light emission 
from STM has been used to reveal plasmonic properties of electrons confined to triangular 
islands on surfaces of noble metals 1231 . All these studies clearly show that, at least at low 
temperatures, quantum confinement plays an important role leading to inhomogeneous electron 
densities at interfaces that differ from the piecewise constant bulk values. 

2 FROM BULK APPROXIMATION TO A LOCAL DENSITY 
APPROXIMATION 

The qualitative optical consequences of an inhomogeneous electron density can be explored 
within Lindhard's "local density approximation" for the dielectric function I24II25I 

e{r,iu) = l- '^^ . n(r) (1) 

where n(r) = n{r)/ no accounts for the density variations relative to the bulk density of elec- 
trons no — mujp/e'^, with Up being the bulk plasma frequency. To mimic the Ohmic damping 
we have included a damping rate 7 and in the optical regime 7 ^ w ll26l . While several 
studies have considered the consequences of local density variations for the surface-plasmon 
dispersion relation, the overall effects have been reported to be minor and for planar surfaces 
the surface-plasmon resonance at w = ujp/\/2 appears unaffected by the detailed density pro- 
file at the metal-air surface |27|. However, it remains an open question to which extent the 
detailed field distribution will be affected by a smoothly varying electron density in contrast to 
the discontinuous drop in the electron density assumed by the bulk approximation. 

Zero-index phenomena in metamaterials are receiving increasing attention ||28H33l , and the 
recent results by Litchinitser et al. fSTl show that for a metamaterial with the effective constitu- 
tive coefficients varying linearly, e(z) cx z and /i(z) cx z, a pronounced field enhancement will 
occur at z = where e = /i = 0. In this paper we show that plasmonic structures employing 
Nature's own metals may support large field enhancement at the zero-e surface appearing in the 
vicinity of the air-metal interface, where the real part of the polarization of the dilute electron 
gas will exactly balance the polarization of the vacuum background. As shown below, the par- 
ticular position of the zero-e surface depends on the frequency oj of the incident radiation, and 
may be either in air or in the metal, depending on the local geometry. 

3 PHENOMENOLOGY OF FIELD ENHANCEMENT 

The basic mechanism behind the field enhancement can be understood from the continuity 
condition for the electrical displacement field, which stipulates that n ■ D ^ e{n ■ E) be 
continuous (here n is the normal vector of a surface element). At the zero-e surface the normal 
component E± = n ■ E thus has to diverge to keep the product eE± finite. 

From Eq. ^ it follows that in the regions where Re(e) vanishes, there will still be a small, 
but finite imaginary part: Im(e)r{o(e)=o = l/'^- For an electrical field Eq incident at an angle 9 
we thus expect that the normal component n-£^o = Eq sin 9 is enhanced by a factor proportional 
to 1/7 on the zero-e surface. The intensity enhancement can then be estimated as 

\E\^/\Eo\^ (xsm^9/-/'^ (2) 

We emphasize that this mechanism of field-enhancement is generic and robust as it only depends 
on the assumed continuous negative-to-positive variation of the dielectric function at the metal- 
air interface, in combination with a Maxwell boundary condition ll35l . rather than on the details 
of the excitation of strongly localized surface-plasmon resonances with pronounced sample-to- 
sample fluctuations. 



4 EXAMPLE 1: JELLIUM MODEL FOR THE METAL SURFACE 

The low-temperature density profile for Jellium is a function of rg/oo only, with the radius 
of the free-electron sphere = (3/477710)^^"^ and ao ~ h? /me^ being the Bohr radius. The 
position-dependent electron density may be obtained self-consistently with the aid of density- 
functional theory |f36l . thus taking into account both exchange and correlation effects in the 
inhomogeneous electron gas. In Ref. ||36| , results are tabulated for a range of values of rg/aQ. 
The resulting density n{z) is shown in Fig.[Tla) for rg/ao = 2, a typical value for gold or silver. 
As seen, n(z) = n{z)/no varies continuously on the length scale of the Fermi wavelength 
rather than exhibiting a discontinuous drop at the surface of the metal (z = 0). In the metal 
{z > 0), there are Friedel oscillations near the surface, while for z < the decay of the density 
is caused by the finite value of the work function W, which allows electron wave functions to 
penetrate into the classically forbidden region. According to Eq. ([T]i, the dielectric function will 
inherit the Friedel oscillations in the metal, as well as approach the value in vacuum outside the 
metal surface. As a consequence, the dielectric function passes through zero just outside the 
metal surface [green line in Fig.[TJa)]. Similar spatial variations of the dielectric function have 
been reported in more elaborate studies of e.g. semiconductor quantum dots j;37| and recently 
also in doped semiconductors ||38l . 

In the case of A = 2nc/uj = 633 nm the zero-e surface appears at a distance of 0.25 nm 
from the surface (z ~ 0). Fig. [Ttb) illustrates the field enhancement in the case of a plane 
wave incident at an angle of 9 ~ 60°. The numerical results were obtained by solving the wave 
equation 

^2 

W xW X E(r) = e(r,u) — E(r) (3) 

with the aid of a finite-element method employing an adaptive mesh algorithm (Comsol Multi- 
Physics). A comparison with e(z) in Fig.[T](a) shows that the field enhancement indeed occurs 
at the point where e = 0. We emphasize that this phenomenon has no counterpart within the 
classical bulk approximation of the dielectric function. For a realistic value of the damping 
(7 = 1 X lO^^rad/s) the intensity may be enhanced by more than 3 orders of magnitude, 
which suggests that in a SERS experiment the Raman signal could be enhanced by 6 orders 
of magnitude |39| without the need for any geometry-induced localized resonances. To further 
test the predictions of Eq. (|2|i we have explored the dependence on the angle of incidence as 
well as the particular value of the plasmon damping. Fig.|2](a) illustrates how the peak height 
depends on damping. The enhancement is limited by the plasmonic damping with the slope of 
the dashed lines being in full accordance with the power (-2) of 7 in Eq. Panel (b) illustrates 
the angle dependence with the dashed line showing the sin^ 9 dependence of Eq. (O, indicating 
a pronounced enhancement for oblique incidence. 

5 EXAMPLE 2: A METALLIC NANOSTRUCTURE 

We next explore the field enhancement in the presence of strong confinement of electrons in 
metallic nanostructures. Here we limit ourselves to a model of non-interacting electrons in wire 
geometries with nanoscale cross sections so that the problem effectively is two-dimensional. 
The energy states of the electron system are then given by S'{k) = + (f^ with the trans- 

verse energy components being quantized and governed by a two-dimensional Schrodinger 
equation 



2m 



yly + V{x,y) 



i'^{x,y) = S^%l^^{x,y). (4) 



For the confinement potential V we adjust the height in accordance with the work function of 
the metal W = Vq — /u, where /u is the chemical potential. In equilibrium the states are populated 



according to the Fermi-Dirac distribution function and we find 



n{x, y) = = ^ 'y- r (5) 



no 



where A = J^dxdy is area and P{x) = Lii/2[^ cxp(— a;)] with Li„(2:) being Jonquiere's 
polylogarithm function. Equilateral triangular nano islands have been the subject of numerous 
investigations I140II41I and in Fig.[3]we consider field enhancement in such a structure. Fig.[3ja) 
shows results within the bulk approximation, where an artificial rounding of the corners has 
been added to circumvent the effect of an otherwise divergent electrical field 121. Clearly, field 
enhancement is modest since no pronounced surface plasmon resonances are excited. Fig. [3b) 
on the other hand is based on the spatially continuously varying dielectric function of Eqs. ([T]i 
and ^ and shows a pronounced field enhancement along the zero-e surface. We also note that 
while the geometry has arbitrarily sharp corners, the electron density and dielectric function 
themselves are still smooth, showing that quantum confinement provides a built-in smoothing 
even in the presence of sharp geometrical features. Fig. Ua) and (b) show in more detail the 
field enhancement along the lines indicated in Fig. Ob). For the results in Fig. |4|b), the peak 
height clearly matches the one-dimensional results in panel (b) of Fig.|2]for an incident angle 
of 6* = 30°. Similarly, the result within the bulk approximation resembles typical results, see 
e.g. Ref. 135)1 . i.e. the field enhancement is more modest. In Fig.|4{a), the incident wave propa- 
gates almost parallel to the zero-e surface [6 ~ 90°), thus causing a peak that is even higher at 
the two right-most corners of the triangle. On the other hand, at the left-most corner the wave 
has normal impact on the zero-e surface (9 ~ 0°) so that no pronounced field enhancement is 
observed. It is interesting to note that the enhancement may occur either outside [as in Fig. [Tib)] 
or inside [as in Fig.[3lb)] the metal, depending on the local geometry and the frequency relative 
to the plasma frequency. 

6 CONCLUSIONS 

We predict that a large non-resonant field enhancement may occur at zero-epsilon surfaces near 
the air-metal interface. Theoretically, the enhancement occurs in a wide class of models where 
the discontinuous jump of the dielectric function of the usual bulk description of at an air-metal 
interface is replaced by a continuous function e(r, w) that has the bulk values as limiting values 
away from the interface. We discussed two microscopic models for such a continuous variation 
at the interface. The predicted field enhancement has therefore no counterpart in the classical 
bulk treatment of plasmonic field enhancement, though it is intimately related to the familiar 
boundary condition for the electrical displacement field. In particular, for hght impinging on an 
interface we discussed the dependence of the field enhancement on the angle of incoming light 
and on the damping in the metal. Furthermore, the enhancement constitutes another mechanism 
for the enhancement of Raman signals at metal interfaces (SERS). 

Interestingly, Feigenbaum et al. have recently reported their experimental observation of 
zero-index driven field enhancement at the surface of conducting oxides (i.e. with lower electron 
density compared to good metals) |[38l . We emphasize that the decay of the electron density 
into vacuum (or more generally into a dielectric) could have important implications also for 
nano-gap structures (such as the gap formed between the corners of two nearby triangles BOl ). 
where quantum tunneUng of electrons will introduce a small but finite electron density inside the 
gap, thus also modifying the local dielectric function and the ability to support short-circuiting 
currents 1191. 
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Fig. 1 . (a) The blue line shows the position dependent electron density within the Jellium model 
(rs/ao = 2). The green line shows the corresponding derived dielectric function, (b) The field 
enhancement for a plane wave incident on the surface at an angle of = 60°, for A = 633 nm, 
= 1.37 X 10^*' rad/s, and 7 = 1.0 x 10^"^ rad/s. 
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Fig. 2. (a) Maximal field enhancement as a function of damping for parameters of Fig. [lib). 
The dashed lines indicate a 7^^ dependence, (b) The peak height as a function of incident angle 
for 7 = 1.0 X 10^^ rad/s. The dashed line indicates a sin'^ 9 dependence, Eq.©. 




Fig. 3. (a) Intensity enhancement (log scale) for a wave incident from the left on a metallic 
equilateral triangle, computed with the bulk dielectric function, but with rounded corners (radius 
of 5 nm). (b) Field enhancement for the inhomogeneous density. Field profiles along the lines 
ri —7- r2 and — )• are shown in Fig.|4] 
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Fig. 4. Field enhancement along the lines defined in Fig. 3. (a) ri — > r2- (b) — > r^. Note that 
the integrated density is the same for the two structures in Fig. [3] thus causing slightly different 
densities and field suppression inside the metal structures (compare dashed and solid lines for 
the inhomogeneous and homogeneous cases, respectively). 



